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FOREWORD 


The  work  presented  in  this  report  was  performed  by  Kaman 
Aerospace  Corporation  under  Contract  DAAJ02-70-C-0012 
(Task  1F162204AA4301)  for  the  Eustis  Directorate,  U.  S.  Army 
Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Virginia.  The  program  was  implemented  under  the 
technical  direction  of  Mr.  Joseph  H.  McGarvey  of  the 
Reliability  and  Maintainability  Division* and  Mr.  Arthur  J. 
Gustafson  of  the  Structures  Division.’MrlThe  report  is  pre¬ 
sented  in  four  volumes,  each  describing  a  separate  phase  of 
the  ba  ,ic  theory  of  structural  dynamic  testing  using  impedance 
techniques . 

Volume  I  presents  the  results  of  an  analytical  and  numerical 
investigation  of  the  practicality  of  system  identification 
using  fewer  measurement  points  than  there  are  degrees  of 
freedom.  The  parameters  in  Lagrange's  equations  of  motion, 
mass,  stiffness,  and  damping  for  a  mathematical  model  having 
fewer  degrees  of  freedom  than  the  linear  elastic  structure  it 
represents  may  be  determined  directly  from  measured  mobility 
data.  Volume  II  describes  the  method  of  system  identification 
wherein  the  necessary  impedance  data  are  experimentally  deter¬ 
mined  by  applying  a  force  excitation  at  a  single  point  on 
the  structure. f  Volume  III  presents  a  method  of  determining 
the  free-body  dynamic  responses  from  data  obtained  on  a  con¬ 
strained  structure.  Volume  IV  describes  a  method  of 
obtaining  the  equations  for  the  combination  of  measured 
mobility  matrices  of  a  helicopter  and  its  subsystems.  The 
response  of  the  combination  of  a  helicopter  and  its  sub¬ 
systems  is  determined  from  data  based  on  the  experimental 
results  of  the  main  system  and  subsystems  separately. 


^Division  name  changed  to  Military  Operations  Technology  Division. 
^♦Division  name  changed  to  Technology  Applications  Division. 
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INTRODUCTION 


The  success  of  a  helicopter  structural  design  is  highly 
dependent  on  the  ability  to  predict  and  control  the  dynamic 
response  of  the  fuselage  and  mechanical  components.  Con¬ 
ventionally,  this  involves  the  formulation  of  intuitively 
based  equations  of  motion.  Ideally,  this  process  would 
reduce  the  physical  structure  to  an  analytical  mathematical 
model  which  would  predict  accurately  the  dynamic  response 
characteristics  of  the  actual  structure.  Obviously,  the 
creation  of  such  an  intuitive  abstraction  of  a  complicated 
real  structure  requires  considerable  expertise  and  inherently 
includes  a  high  degree  of  uncertainty.  Structural  dynamic 
testing  is  required  to  substantiate  the  analytical  results. 
The  analysis  is  modified  until  successful  correlation  is 
obtained  between  the  analytical  predictions  and  the  test 
results . 

This  report  describes  the  theory  of  structural  dynamic 
testing  using  impedance  techniques  as  applied  to  a  mathe¬ 
matical  model  having  fewer  degrees  of  freedom  than  the 
structure  it  represents.  The  test  information  is  obtained 
with  single  point  excitation  of  the  model.  Reference  1 
describes  the  method  of  obtaining  a  model  directly  from  test 
measurements  for  a  hypothetical  structure  which  has  the  same 
number  of  degrees  of  freedom  as  the  mathematical  model.  In 
reality,  the  number  of  degrees  of  freedom  of  a  physical 
structure  is  infinite/  therefore,  the  usefulness  of  model 
identification,  necessarily  with  a  finite  number  of  degrees 
of  freedom,  using  impedance  testing  techniques  depends  on 
the  ability  to  simulate  the  real  structure  with  a  small 
mathematical  model.  Reference  2  illustrates  the  method  of 
obtaining  a  model,  using  impedance  testing  techniques,  that 
is  comprised  of  less  degrees  of  freedom  than  the  physical 
structure  it  approximates.  That  method  required  measured 
mobility  data  obtained  at  selected  points  of  the  structure 
with  the  force  input  applied  at  each  of  the  prescribed 
locations.  The  present  theory  is  similar  to  that  of  Ref¬ 
erence  2  except  that  the  excitation  is  applied  at  only 
one  point  on  the  model,  thereby  substantially  reducing  the 
mobility  data  essential  to  the  analysis. 

The  process  of  deriving  the  equations  of  motion  from  test 
data  is  referred  to  as  system  identification.  The  only 
input  information  required  in  this  theory  is  measured 
mobility  data  obtained  with  the  excitation  at  only  one  point 
on  the  model  and  the  approximate  natural  frequency  of  each 


mode.  This  information  can  be  readily  obtained  from  impedance 
testing  of  the  actual  structure  over  the  frequency  range  of 
interest  yielding  the  second  order,  structurally  damped 
linear  equations  of  motion. 

System  identification  theories  of  any  practical  engineering 
significance  must  be  functional  with  a  reasonable  degree  of 
experimental  error.  In  this  report,  a  series  of  computer 
experiments  incorporating  experimental  errors  was  documented. 
This  report  presents  an  extension  of  the  analysis  derived  in 
Reference  2  whereby  an  identified  model  with  a  finite  number 
of  degrees  of  freedom,  obtained  from  impedance  type  testing 
with  excitation  at  only  one  point  on  the  structure,  simulates 
the  actual  structure  wherein  the  number  of  degrees  of  freedom 
is  infinite. 
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DERIVATION  OF  THE  SINGLE-POINT  ITERATION  PROCESS 


As  indicated  in  References  1  and  2,  the  mobility  of  a  struc¬ 
ture  is  given  by 


l*j  -  •*!  fYi  (tol-l  1*1 


0) 


(1) 


K 


K 


With  excitation  at  station  k,  the  responses  at  station  j, 
including  k,  are  obtained.  These  provide  the  k-th  column 
of  the  mobility  at  a  particular  forcing  frequency  : 


N 


{Yj001}  =  Yi<l)*ki{*}i  - 


(2) 


l<j<J,  l<i<N 


This  represents  a  column  of  mobility  values,  each  element  of 
viiich  is  the  response  at  a  point  of  interest  on  the  structure 
with  excitation  at  station  k  and  at  forcing  frequency 

Similarly,  with  the  exciter  remaining  at  station  k,  the  k-th 
column  of  the  mobility  at  another  frequency,  w2  ,  can  be 
obtained  : 


N 


{Yj(k)2>  *  Yi<2)*ki<*}i  *  '®>{ViAi} 


(3) 


The  mobility  columns  represented  by  (2)  and  (3)  may  be 
combined  into  one  matrix: 


[{Yj(k)lHYj(k)2  }]  -  1*1  CCY u^iHY-Ai}] 


Jx2 


"  l*]  r<|,kiJ  I{Y±1}{Yi2>] 


JxN  NxN 


Nx2 


In  general,  for  P  forcing  frequencies  (1  <_  p  <_  P) , 
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JxP 


JxN  NxN  NxP 


(4) 


(5) 


If  J  >  P,  Equation  (5)  is  a  set  of  more  equations  than  un¬ 
knowns  for  which  there  is  no  solution.  Equation  (5)  can  then 
be  written  as 


j  OOP1 

- 

+  [R. 
IP 

JxP 

JxN  NxN  NxP 

Jxp 

(6) 


where  R^p  is  the  residual  associated  with  the  j-th  station 
and  theJp-th  forcing  frequency. 

As  described  in  References  1  and  2,  the  imaginary  displace¬ 
ment  mobility  contains  significant  information  relating  to 
modes  associated  with  natural  frequencies  in  proximity  to 
the  forcing  frequency.  As  shown  in  Reference  3,  accurate 
estimates  of  the  modal  vectors  may  be  obtained  by  considering 
only  the  effects  of  modes  proximate  to  the  forcing  frequency. 
Therefore  the  analysis  will  employ  only  Q  modes,  where  Q  is 
less  than  N.  Consider  the  imaginary  displacement  mobility 


[Yj(k)p]  “ 


[4]  [Y*1]  +  [R.  1 

i  lTki  ip  ]p 


(7) 


TY  •  1 

The  dominant  element  in  each  row  of  the  1  ip  matrix  will  be 
the  modal  mobility  measured  at  the  forcing  frequency  in 
proximity  to  a  particular  natural  frequency.  Normalizing 
the  rows  of  the  aforementioned  matrix  on  the  largest  element 
yields 


‘V  - 


(8) 


*1 

Y  . 

where  in  is  the  maximum  value  of  the  l-th  row. 
may  be  rewritten,  incorporating  Equation  (8)  : 


Equation  (7) 


[Y*  .  1  =  [4]  n>i,-Y*I4[S.  ]  +  [R.  ] 

](k)p  Tki  in  ip  3P 


(9) 


The  [Sip]  matrix  can  be  evaluated  by  considering  the  ex¬ 
pression  for  the  imaginary  displacement  modal  mobility 
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Therefore  from  Equation  (8), 


S. 

ip 


2  2 

o  0) 

Si  +  (1 - %) 

ni 

2  a)  2  2 

gi  +  (1  -  -J2y> 

n. 


(10) 
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Because  g^ ,  the  structural  damping  coefficient  of  the  i-th 
mode,  is  generally  quite  small,  typically  of  the  order  5 
percent,  the  [S]  matrix  can  be  accurately  estimated  by 
assuming  g.^  =  0,  thus,  requiring  knowledge  of  only  the 
forcing  frequencies  and  the  natural  frequencies.  It  will  be 
shown  that  an  accurate  estimate  of  S  is  not  necessary, 
although  helpful,  as  iterations  will  converge  on  the  best 
values  in  S  in  the  least-squares  sense. 


The  matrix  Equation  (9)  has  no  solution.  An  approximation 
to  a  solution  may  be  defined  as  that  which  makes  the 
Euclidian  norm  of  the  matrix  of  residuals  a  minimum.  This, 
as  will  be  proved  later,  is  given  through  use  of  the  pseudo¬ 
inverse. 


Equation  (9)  will  be  solved  utilizing  matrix  iteration 
techniques  using  ,<.(0),  as  a  first  estimate.  As  indicated 

lbip  J 

in  the  following  sections,  the  modal  vector  matrix  with 
respect  to  which  the  Euclidian  norm  of  the  residuals  is  a 
minimum  is  given  by 


q)kiYin 

where  [s|p^]+  is  defined  as  the  generalized  inverse  or 
pseudoinvirse  of  [S  ^  ]  and  is  given  by 
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(13) 


where 


[s.(0)]  =  [s.(0)  ]T( [s.(0)  ]  [s.(0)  ]T) 

ip  1  1  ip  J  ' 1  xp  J  ip 


fs  £>ns£)]  =  uL] 


It  follows  then  that 


[Y1,,.  ]  =  ]  IV. -Y?  *4  [S.(0)  ]  +  [R  ] 

j(k)pJ  J  lYki  in  ip  JP 


(14) 


in  which  the  Euclidian  norm  of  [rP^]  is  a  minimum  with 
respect  to 

Using  [^^],a  matrix  [S^p^  ]  can  be  found  to  give  an  equation 


tY?(k)pi  -  tt(1)mkiY*;jts^,i  +  [R<1>] 

such  that  the  Euclidian  norm  of  [Rj^  ]  is  a  minimum  with 
respect  to  [S  ] .  This  is  given  by 


(15) 


[g  (1)  l  -  r _ 1 

lSip  J  r.  „*I 


<J>,  .  Y . 
ki  in 


*rJ[$(1)]  +  [yJ 


j  OOP 


] 


where 


(16) 


m+  =  ( [<t|]T[(|>]  )”1  [<J>] T  and  t<j)]+[<f>]  =  [IR] 


(17) 


It  is  apparent  from  the  first  cycle  of  the  iteration,  by 
comparing  Equations  (11)  and  (15) ,  that  the  process  consists 
of  alternately  dealing  with  the  left  and  right  identity 
matrices.  At  each  successive  iteration,  a  solution  is  found 
that  minimizes  the  Euclidian  norm  of  the  residual  matrix  with 
respect  to  the  newly  found  matrix  of  either  [S]  or  [$] . 
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'«BWW»^W» 


In  simplified  notation,  the  q-th  iteration  becomes 
[*(q)]  =  [Y1]  [S(q~1)  ]•*>■■  1  Vf  J 

♦kiYii; 

and 

[S(q)]  =  ^  1 . TjJ  [»(q)]'*‘[YIl 

kiYin 

The  next  iteration  is 


[$(q+1) ] 
[S(q+1)]  = 


=  yI[S(q)]+r— 

*kiYin 

r — nfYj[4(q+1)]+fYI] 

^iYi; 


(18) 


(19) 


This  is  the  basic  algorithm  used  in  the  matrix  iteration 
procedure. 


DETERMINING  THE  MODAL  PARAMETERS 


From  Equation  (6)  of  the  previous  section,  one  column,  which 
is  at  a  particular  forcing  frequency,  p,  with  the  excitation 
at  station  k,  can  be  written  as 


fYj  <kp)  1 


'W 


(20) 


The  number  of  modes,  Q,  included  in  Equation  (20)  cannot 
be  greater  than  the  number  of  points  of  interest  on  the 
specimen,  J,  and  generally  will  be  much  less  since  only  those 
modes  which  have  significant  effect  on  the  mobility  at  the 
forcing  frequency,  u)p,  wij.1  be  considered.  Ordinarily,  the 
number  of  modes  used^will  not  be  greater  than  3  or  4  for  any 
given  forcing  frequency,  and  these  will  be  the  modes  in  the 
vicinity  of  the  forcing  frequency  in  question. 

The  real  and  imaginary  modal  mobilities  are  calculated  from 


and 


(Y 


*R 
i  (p) 


}  =  r-r^-J  [*]+{Y? 


^ki 


j  (kp) 


{Yi(P>}=  rf[*1  +  {YW 


(21) 

(22) 


From  Reference  1  the  real  displacement  mobility  can  be 
calculated  as 


*n  i  i-o)  2/n.2 

*R  _  1  _  p  '  l 

iwP  g.2  +  (1  -  wp2/Di2)2 


(23) 


and  the  imaginary  modal  mobility  by 


, 

,,  2  /n  2.  2 

(1  -  u,p  /fl.  ) 


The  real  modal  impedance  can  be  written  as 


(24) 
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(25) 


,*R 


Y*R 

1U)P 

- *R - 2 - - 3Ff - 7 

(y.k  r  +  (y.1  r 

lwp  la)p 


Substituting  Equations  (23)  and  (24)  into  (25)  yields 


Z*R  =  K-  (1  "  w  2/n.2)  (26) 

10)  '  l  p  '  1 

P 

From  Equation  (26)  it  is  observed  that  the  modal  impedance 
is  a  linear  function  of  the  square  of  the  forcing  frequency. 

The  forcing  frequency  at  which  the  modal  impedance  becomes 
zero  is,  therefore,  the  natural  frequency.  From  a  least- 
squares  analysis  of  modal  impedance  as  a  function  of  forcing 
frequency  squared,  proximate  to  the  natural  frequency,  the 
generalized  stiffness  of  the  i-th  mode  and  the  natural  fre¬ 
quency  of  the  i-th  mode  can  be  calculated. 

The  generalized  mass  associated  with  the  i-th  mode  is  given 
by 

Hi  =  Vfii2  (27) 


The  structural  damping  coefficient  may  be  determined  from 

2  Y*1 

gi  =  (-£2  _1)  ~"*R 

V  Yiwp 


(28) 


EQUATIONS  OF  MOTION 


There  are  two  basic  types  of  dynamic  mathematical  models 
describing  structures.  The  conventional  type,  covering  as 
many  modes  as  degrees  of  freedom,  is  called  "Complete  Models" 
and  is  considered  in  References  1  and  2.  The  other  type 
labelled  "Incomplete  Models"  considers  fewer  modes  than 
points  of  interest  on  the  structure  and  was  first  described 
in  Reference  5.  Using  the  methods  described  herein,  it  is 
possible  to  identify  either  a  complete  model  or  a  form  of 
incomplete  model. 


Incomplete  Models 


Consider  a  rectangular  identified  modal  matrix  which  has  J 
rows  indicating  the  points  of  interest  on  the  structure  and 
Q  columns  representing  the  modes  being  considered  where  J  >  Q. 
The  influence  coefficient  matrix  for  the  incomplete  model  is 
given  by 

[Cinc]  =  Tjq-J  [4>]T  (29) 

The  above  matrix,  similar  to  all  incomplete  model  parameter 
matrices,  is  singular,  being  of  rank  Q  and  order  J.  The 
mass,  stiffness  and  damping  matrices  for  the  incomplete 
model  are 


tm.  i  =  m+TVj  [*]+ 

xnc  i 

[k.  j  -  m+TrK.4  m+  oo) 

me  l 

Mine1  -  m+TrgiKiJm+ 

The  classical  modal  eigenvalue  equation  has  the  analogous 
incomplete  form 


i 


Complete  Models 


(31) 


For  the  complete  model  the  identified  modal  vector  matrix 
is  square,  having  the  same  number  of  degrees  of  freedom  as 
mode  shapes;  that  is,  J  =  Q.  The  influence  coefficient  matrix 
is  given  by 
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vmrrmmmmm 


7*#*nmKm *n9m*. 


[c]  = 


m  ti/iy  mA  = 


N  1  T 

i=l  *i  1  1 


(32) 


The  mass,  stiffness  and  damping  matrices  for  the  complete 
model  are 


[m]  =  [<j>] 

[k]  =  [4>]‘TfKiJ  W"1 
[d]  =  t<l>3  “T  tgi  Ki'l  1*1 


as  indicated  in  Reference  1. 

Full  Mobility  Matrix 

The  full  mobility  matrix  of  either  complete  or  incomplete 
models  is  given  by 


[Y]  =  [$]fY*J[$]T  (34) 

where  for  the  complete  model  the  [$]  matrix  is  square,  having 
J  columns  and  J  rows.  However,  in  the  case  of  the  incomplete 
model  the  modal  matrix  [$]  is  rectangular,  having  J  rows  and 
Q  columns,  where  J  >  Q. 
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PROOF  THAT  THE  PSEUDOINVERSE  MINIMIZES  THE  NORM 
OF  THE  RESIDUALS 


Take  the  transpose  of  Equation  (9)  and  write  the  equation 
for  one  column  of  the  transpose  of  the  mobility  matrix: 


IYWT=  <vt<vt  +  [vT 


{Y 


(j)Op 


>  «  [SipiT{*(j)i}  +  (rm^> 


(j)p 


(35) 


{R(j)p} 


(Y 


(jk)p}  "  [Sip]  {<t,(j)i} 


(36) 


{R(j)p}  {R(j)p}  *  {Y 


(jk)p}T{Y(jk)p}  "  {Y 


(jk)p}TlSip]TU(j)i}  ' 


lSip^Y(jk)p^  +  IS,„][S4J 


(j)iJ  l°ipJli3ipJ  M  j)  i 


(37) 


Equation  (37)  is,  of  course,  a  scalar  product  and  it  is 
recognized  that  the  derivative  of  a  scalar  with  respect  to 
a  vector  is  a  vector;  in  other  words,  Equation  (36)  is  a 
vector  in  p-dimensional  space  and  Equation  (37)  is  its  dot 
produce  on  itself  -  that  is,  its  length  squared.  We  wish  to 
find  the  vector  { }  which  makes  the  length  of  the  residuals 
vector  a  minimum. 


Take  the  partial  derivative  of  Equation  (37)  with  respect  to 
T 

{ ^  ( j )  i )  an<*  set  e(3ual  to  zero  to  obtain  the  modal  vector 
for  which  the  Euclidian  norm  of  the  residuals  is  a  minimum: 


0  -  -2[S^>HVtjk)p)  + 


} 


or 


U(j)i}  (iSip  “ip  3  3  Sip  ,1Y(jk)p' 


(38) 


and 


rad)  |y  (0)  i  T . (0)  (0)  .  Tv -1 

{*(j)i}  tY(jk)p3  lSip  3  Sip  Sip  J  3 


(39) 
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as  the  inverted  matrix  is  symmetrical.  Equation  (39)  is  any 
row  in  Equation  (12).  The  sum  of  the  minimum  Euclidian  norms 
of  the  rows  of  a  matrix  is,  by  definition,  the  minimum 
Euclidian  norm  of  the  matrix,  and  it  therefore  follows  from 
Equation  (39)  that 

|t(1)>  *  'Yj(k)p1'sip)>T"sip>HS^>)T)-1 

which  is  given  by  Equations  (12)  and  (13).  Q.E.D.  The  basic 
observation  which  makes  the  above  proof  of  the  pseudoinverse 
possible  should  be  credited  to  Klosterman,  Reference  (4). 

To  show  that  the  [S]  matrix  obtained  using  the  pseudoinverse 
of  [4>]  minimizes  the  norm  of  the  residual,  write  the  equation 
for  a  column  of  Equation  (9)  : 


{Yj(kp)}  “  [#1{Si(p)}  +  {Rj(p)} 


{Rj(p)  }  *  {YjOcp)}  "  [*1{Si(p)} 


(40, 


{Rj(p)}T{Rj(p)}  "  {Yj(kp)}T{Yj(kp)}  "  {Yj(kp)}Tm{si(p)} 

-  {si(p)}TmT{Yj(kp)}  +  {si(p))TmTm{si(p)} 


(41) 


3{Rj #n\ )  (Rj /n\ }  (i )  . 

Set  - li£2_ —  JAP?  -  -  o  and  solve  for 


91si (P»> 


(sMp)>  -  t»lT<*Jtkp)> 


(42) 


or 


ts^’i  -  n*iTi*J>‘1it]TirJ(lt)p) 


(43) 


which  is  the  same  as  Equation  (16).  Q.E.D. 
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PROOF  THAT  ITERATIONS  USING  THE  PSEUDOINVERSE  OF  S  AND  $ 
CONVERGE  MONOTON ICALLY  ON  MINIMUM  SUM  OF  RESIDUAL  SQUARES 


In  the  q-th  iteration,  where  q  is  odd. 


[Y 


jOOp 


3  -  [*(q‘1))[S4(r1)]  +  CRir^l 


’ip 


jP 


(44) 


[*<q)]  =  [Yj(k)p][s|q“1)]+  -  [*(q_1)]  +  [Rj(q“1)]  tS(q_1)]  + 


(45) 


because  [S3  [S]  -  [1^1*  Then 


txj  OOP1  “  t*(q)1<-1>]  ♦  l*g»l  (46) 


Substitute  Equation  (45)  into  Equation  (46): 


"W  -  (4<q"1)iisip"l)i  4- 

or  +  [RiS,)  (47) 

-  i*w  -  <i,> +  <l,)isirl,i+'sr)> 

+  [R  fq)  ] 

Therefore  DP 

[r]^]  -  [Rjp_1)  1  (>H  -  fsip’’1)l+lsip"1)l)  (48) 


JxP 


JxP 


PxP 


The  p-th  row  of  [R.jq^ 1  is 


fR(q)  1T  -  fR(q-1)Vr(>U 
(Rj(p)}  -  (Rj(p)  )  (TU 


Isip‘1>l  +  [si5“1)l) 


,R1(<(P)  ,T(r3(<(p)  1  -  tRjqp)))T(>Ii  -  'sip'1>1+,si?'1)  1,(111 


-  [S 


(q-1)  j  +  [g  (q-1) ,  ^T/o  (q-1) 


iP 


V4~x/  i.irR  -x/  > 

ip  J)  1Rj(p)  * 
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But  [I]  -  [S.(q'1)]  +  tS.(2"1)  is  symmetrical  and,  from  Equation 
(13),  ip  ip 


ts{q"1)](s{q"1)3+  -  IIL1.  Therefore, 


lRj(p)}  tRj(p)}  lRj(p)  }  lRj(p)  j 


_  (RMiT.Jq-D.  +  fJq-DwjCq-l), 

tRj(p)  }  [Sip  ]  Sip  ]iRj(p)  } 


(49) 


[S^q  maximally  ranked  in  its  rows,  of  rank  Q  where 

1  _<  i  _<  Q.  Therefore  [sfq  ^  ]  lS^q  ^]T  and  its  square  root 

( tS±p  ^Sip  ^]T)1//2  are  nonsingular  of  rank  Q  and 

symmetrical.  Now,  [S.^q  ^]  +  [S^q  is  real,  symmetric  and 

ip  rp 

singular.  It  is  known  that  a  real  symmetric  matrix  [A]  of 
rank  Q  is  positive  semidefinite  if  and  only  if  there  exists  a 

matrix  [C]  of  rank  Q  such  that  [A]  =  [C]  [C] .  Let 

( ]T)-1/2  [sJq-1^  =  [C] ,  rectangular  of  rank  Q. 


ISip’1)lTnS^“1)][S^"1)]T)“  2  ([s{q"1))[s{q"1)]T)  I[s{q'1)] 


CTC 


lsip  1>  jT( Isip”1> 1  Isip-1)]T)_1 lsip’1>  1 

-  Csg^Vcsg-U] 


(50) 


Therefore  [S^q  ^^]+(s|q  is  positive  semidefinite  and 

}T [s 5q_1  ^  ]+  [sJq-1^  ]  {r!c|"^  }  in  Equation  (49)  must  be  a 
3  iP  rp  3  »P 

nonnegative  number.  But  the  first  term  on  the  right  side 
and  the  left  side  of  Equation  (49)  are  also  necessarily 
nonnegative.  Therefore 


(51) 


Por  the  alternate  calculation,  q  odd 

Is(q>,  -  [♦''J>]+i^(k)p] 
But  tY;(10p]  -  s£-l>i  +  i «j?i'  80 

ts‘|> ]  -  +  I*<,>)+|Rj(5>] 

Substituting  IS*’1)  for  |s|^-1> ) .  we  obtsin 


(18) 

(46) 

(52) 


[* 


(q)ns  <2>]  ♦ 1 


•ip  *  ■  i*'jP 

-  t*(q))is{q'1)]  +  [*(q))(4(q)]  +  tRjp)l  +  tRjp+1)l 

(53) 

Prom  Equations  (46)  and  (53)  , 

-  «*W  -  «8?>  +  +  <"> 


or 

[R'q+1)]  -  OH  - 
DP 

Compare  Equation  (54)  to  Equation  (48)* 


(54) 
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Consider  a  column  of  Equation  (54)  {R.^+^}.  Because  of 
Equation  (18),  3  (p) 

8/R(q+i)  >TfR(q+D> 

__j  (p)  J  T  (p)  }  _  A 


*i  (p) 


{Rj1p))}T{Rj^)}  "  {Rjitp)}T(m  '  t*(q)n*(q)]V([i) 

19  Jl*  J  H  Rj(p)'  tRj(p)}  {Rj(p)} 


because  [$]+[$]  =  [IR]  (Equation  15)  and  [4^][$^]  +  is 

T 

symmetrical.  Now  [$  ^  ]  [4>  ^  ]+  =  [$ ^  ]  ( [$  ^  ] )  ”  J 

([$(q^j)  1/2[$^)]T  and  [$^^]  is  necessarily  maximally 

column  ranked.  Therefore,  I  $  ^  ]  [4>  ^  is  positive  semi- 
definite.  The  left  side  of  Equation  (55)  is  the  positive 
difference  between  two  positive  numbers,  and  it  follows 
that 


I  E  (p!‘j+1>>2  <  E  E  tR-ln)>2 

j-1  p-1  3P  j-1  p-1  JP 

l~>o ) 

Equation  (51)  shows  that  the  Euclidian  norm  of  residuals  with 
odd  index  q  is  less  than  the  norm  of  residuals  of  index  q-1; 
Equation  (56)  shows  that  the  norm  of  residuals  of  index  q+1 
is  less  than  the  norm  of  residuals  of  index  q.  Equations  (51) 
and  (56)  show  that  it  is  immaterial  whether  q  is  odd  or  even. 


e  i  (aiq+1))2  <  z  i  (R.(q))2  <  Z  £  (Rfi-”)2 
j-1  p-1  JP  j-1  p-l  jP  j-i  p-i  jP 

Equation  (57)  covers  a  complete  iteration  cycle.  Q.E.D. 
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NOTE  ON  THE  DERIVATIVE  OF  A  SCALAR  WITH  RESPECT  TO  A  VECTOR 

Let  [S]  be  a  square  matrix  of  order  R 

R  R 

{x>T(sHy}  -  i  e  si,x1y^ 

i-l  j-1  13  1  3 


R  R 

{y}T[s]T{x>  -  e  E  s.^x., 

i-i  j-i  J1  1  J 


jLk).TlsIl3gI  -  ?  s.  .y .  -  Is] {y } 

3{x>T  j-1  j  j 


-  &{j}{y{S^  ■  iE1  SijXi  -  tSlT(x> 


R  R 

3{y)T[S] (y)  =  3 (v)T CS]T{y }  ,  a  i-l  j-1 
a{y)T  3{y)T  3{y>T 


afrfjglkl  -  IS]{X}  +  tS]T{X}  -  (IS]  +  (S]THx) 
3(x)T 


) 


i 

1 


j 
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IDENTIFIED  GENERALIZED  MASSES 


Typical  generalized  mass  identifications  are  shown  in 
Tables  I  through  VI.  Table  VII  describes  the  various  models 
for  which  data  is  presented  in  Tables  I  through  VI.  Table 
VIII  presents  a  lumped  mass  description  of  the  twenty-point 
specimen  which  was  used  to  generate  the  simulated  experi¬ 
mental  data.  The  model  stations  used  in  the  various  models 
refer  to  the  corresponding  stations  in  the  twenty-point 
specimen.  Table  I  presents  results  for  model  5C,  which  are 
typical  of  the  results  obtained  for  other  five-point  models. 

Data  are  presented  for  conditions  of  zero  experimental  error 
and  for  simulated  experimental  displacement  mobility  data 
recorded  with  a  random  error  of  +5  percent  and  a  bias  error 
of  +5  percent.  For  the  cases  involving  error,  the  random 
displacement  error  was  computed  using  a  uniformly  distributed 
probability  density  function.  This  error  was  applied  to  both 
the  real  and  imaginary  components  of  the  displacement 
mobility  data.  Table  I  presents  the  effects  of  random  number, 
the  seed  used  in  generating  the  random  error.  The  results 
indicate  the  method  is  extremely  insensitive  to  measurement 
errors  as  applied  herein. 

Table  II  shows  results  for  several  different  five-point 
models.  It  is  apparent  that  no  outstanding  differences 
exist  among  the  models  considered.  The  results  for  the 
twenty-point  specimen,  the  simulated  actual  structure,  are 
also  given  in  the  table  for  comparison.  The  generalized 
mass  distribution  associated  with  each  of  the  models  is  in 
excellent  agreement  with  the  twenty-point  results. 

Tables  III  and  IV  present  results  for  the  nine-point  models 
studied.  Again,  the  calculations  of  the  generalized  masses 
for  the  various  nine-point  models  under  consideration  are  in 
agreement  with  the  simulated  structure. 

Tables  V  and  VI  describe  the  results  of  the  computer  experi¬ 
ments  conducted  employing  the  twelve-point  models.  The 
calculations  produced  acceptable  results  except  for  identifi¬ 
cation  of  the  generalized  masses  of  the  10th  and  11th  modes.  The 
generalized  masses  associated  with  these  models  are  extremely 
small  in  comparison  to  the  remaining  modal  generalized  masses. 
Further,  the  mode  shape  of  the  10th  mode  indicates  lack  of 
response  at  all  points  of  interest  on  the  structure  other 
than  the  first  station.  Therefore,  the  effect  of  the  10th 
mode  is  difficult  to  evaluate  in  the  calculation  of  the 
generalized  parameters. 
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TABLE  I.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
5X5  MODEL*  OF  20  X  20  SPECIMEN 


Computer 

Experiment 

Number 

290 

291 

292 

293 

294 

Random  Disp.  Error 

0 

+5% 

±5% 

±5% 

+5% 

0 

Bias  Disp.  Error 

0 

+5% 

+5% 

+5% 

+5% 

0 

Random  Error  Seed 

- 

5 

13 

421 

1094 

- 

Stations 
(In.)  Mode 

Generalized  Masses 
(Lb-Sec2/ln.) 

1 

8.415 

8.560 

8.543 

8.616 

8.470 

8.534 

2 

4.713 

4.544 

4.619 

4.401 

4.175 

4.449 

3 

.503 

.469 

.493 

.471 

.458 

.495 

4 

1.094 

1.000 

1.050 

1.022 

1.033 

1.087 

5 

.631 

.572 

.651 

.644 

.586 

.630 

22 

32 

43 


*  Model  5C 


**  From  20  x  20  Specimen 


TABLE  II.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
5X5  MODEL  OF  20  X  20  SPECIMEN 


Model 

5A 

5B 

5C 

5D 

Computer 

Experiment 

Number 

296 

297 

292 

295 

Random  Disp. 

Error 

+  5% 

+  5% 

+  5% 

+  5% 

0 

Bias  Disp.  Error 

+5% 

+5% 

+  5% 

+  5% 

0 

Random  Error 

Seed 

13 

13 

13 

13 

- 

Mode 

Generalized  Masses 
(Lb-Sec^/in. ) 

1 

8.544 

8.538 

8.543 

8.568 

8.534 

2 

4.506 

4.506 

4.619 

4.610 

4.449 

3 

.494 

.494 

.494 

.493 

.495 

4 

1.048 

1.047 

1.050 

.994 

1.087 

5 

.653 

.653 

.  651 

.629 

.630 

**  From  20  x 

20  Specimen 

TABLE  III.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
9X9  MODEL*  OF  20  X  20  SPECIMEN 


Computer 

Experiment 


Number 

Random  Bias  Error 

Bias  Disp.  Error 

Random  Error  Seed 

298 

0 

0 

299 

+  5% 

+  5% 

5 

300 

+  5% 

+  5% 

13 

301 

+  5% 

+5% 

421 

302 

+5% 

+5% 

1094 

2** 

0 

0 

Station 

(In.) 

Mode 

Generalized  Masses 
(Lb-Sec2/In. ) 

0 

1 

8.419 

9.283 

9.000 

8.307 

8.253 

8.534 

30 

2 

4.591 

4.462 

4. 350 

4.301 

4.189 

4.449 

140 

3 

.504 

.462 

.472 

.467 

.483 

.495 

160 

4 

1.094 

.975 

1.042 

1.053 

1.095 

1.087 

220 

5 

.631 

.659 

.551 

.577 

.610 

.630 

280 

6 

.761 

.717 

.786 

.674 

.646 

.743 

340 

7 

1.213 

1.152 

1.154 

1.208 

1.052 

1.177 

400 

8 

1.439 

1.371 

1.401 

1.322 

1.370 

1.412 

460 

9 

.813 

.713 

.787 

.860 

.719 

.786 

*  Model  9A 

**  From  20  x  20  Specimen 


i 

t 


22 


TABLE 

•  i! 

> 

H 

IDENTIFICATION 
9X9  MODEL  OF 

OF  GENERALIZED  MASSES 
20  X  20  SPECIMEN 

9 

Model 

9A 

9B 

9C 

20  Pt 

Computer 

Experiment 

Number 

300 

303 

304 

1* 

Random  Disp. 

Error 

±5% 

+  5% 

+5% 

0 

|Bias  Disp.  Error 

i 

+5% 

+  5% 

+5% 

0 

Random  Error 

Seed 

13 

13 

13 

- 

Mode 

Generalizes 

(Lb-Sec" 

[  Masses 
•/In.) 

1 

9.000 

9.015 

9.043 

8.534 

2 

4.350 

4.335 

4.513 

4.449 

3 

.472 

.472 

.472 

.495 

4 

1.042 

1.042 

1.138 

1.087 

5 

.551 

.549 

.584 

.630 

6 

.786 

.783 

.723 

.743 

7 

1.154 

1.243 

1.120 

1.177 

8 

1.401 

1.411 

1.396 

1.412 

9 

.787 

.708 

.791 

.786 

*  From  20  x  20  Specimen 


TABLE  V.  IDENTIFICATION 
12  X  12  MODEL* 

OF  GENERALIZED  MASSES, 
OF  20  X  20  SPECIMEN 

Computer 

Experiment 

Number 

305 

306 

312 

307 

308 

1** 

Random  Disp.  Error 

0 

+  5% 

+  5% 

+  5% 

+  5% 

0 

Bias  Disp. 

Error 

0 

+  5% 

+  5% 

+5% 

+5% 

0 

Random  Error  Seed 

- 

5 

13 

421 

1094 

- 

Station 

(In.) 

Mode 

Generalize 

(Lb-Sec 

d  Masses 

2/In. ) 

0 

1 

8.435 

9.234 

8.474 

8.886 

7.846 

8.534 

30 

2 

4.600 

4.217 

4.556 

4.455 

4.183 

4.449 

60 

3 

.504 

.481 

.488 

.476 

.432 

.495 

120 

4 

1.094 

1.030 

1.150 

1.004 

1.059 

1.087 

140 

5 

.631 

.596 

.596 

.595 

.616 

.630 

180 

6 

.761 

.686 

.722 

.757 

.741 

.744 

220 

7 

1.212 

1.142 

1.182 

1.067 

1.218 

1.177 

260 

8 

1.429 

1.299 

1.232 

1.331 

1.290 

1.412 

300 

9 

.813 

.830 

.797 

.805 

.790 

.786 

340 

10 

.169 

.053 

1.203 

.265 

.565 

.043 

400 

11 

.112 

.091 

.093 

.102 

.120 

.172 

460 

12 

1.135 

1.070 

1.177 

.940 

1.085 

1.050 

*  Model  12B 

**  From  20 

x  20  Specimen 
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TABLE  VI. 

IDENTIFICATION  OF  GENERALIZED  MASSES, 

12  X  12  MODEL  OF  20  X  20  SPECIMEN 

Model 

12B 

12F 

12A 

20  Pt 

Computer 

Experiment 

Number 

312 

311 

309 

1* 

Random  Disp. 

Error 

+  5% 

+5% 

+  5% 

0 

Bias  Disp.  Error 

+  5% 

+5% 

+  5% 

0 

Random  Error 

Seed 

13 

13 

13 

Mode 

Generalized 

(Lb/Sec2 

Masses 

/In.) 

1 

8.474 

8.464 

8.518 

8.534 

2 

4.:- 56 

4.510 

4.492 

4.449 

3 

.488 

.487 

.  487 

.495 

4 

1.150 

1.151 

1.103 

1.087 

5 

.596 

.597 

.  595 

.630 

6 

.722 

.724 

.777 

.744 

7 

1.182 

1.113 

1.159 

1.177 

8 

1.232 

1.242 

1.215 

1.412 

9 

.797 

.743 

.789 

.786 

10 

1.203 

1.043 

-.564 

.043 

11 

.093 

.104 

.0103 

.172 

12 

1.177 

1.119 

1.147 

1.050 

*  From  20  x  20  Specimen 


TABLE  VII.  MODEL  DESCRIPTION 


<N  X  X 


XXX 
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X  XX 


X  X 


VO  X  X  X  X 
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XXX 


X  X 


X  X 


<  £0  U 
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XXX 


XXX 


XXX 


XXX 


XXX 


<  ffl  Cn 
(N  in  (N 


26 


TABLE  VIII.  20-PQINT  SPECIMEN  DES 
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Computer  experiment  309  yielded  a  negative  10th  generalized 
mass.  All  computer  experiments  that  failed  in  this  respect 
gave  drastically  unrealistic  values  of  generalized  mass. 
Ordinarily,  using  different  stations  or  forcing  frequencies 
produced  proper  identification  of  all  modes. 


RESPONSE  FROM  IDENTIFIED  MODEL 


Figures  1  through  12  portray  typical  real  and  imaginary 
acceleration  mobility  response  obtained  from  the  various 
models  considered  in  the  present  study.  In  each  instance, 
the  exact  curve  represents  the  simulated  experimental  data 
for  the  twenty-point  structure,  obtained  with  zero  error. 
Figures  1  and  2  provide  the  effect  of  random  number  seed 
for  a  typical  five-point  model.  Figures  3  and  4  present 
the  results  obtained  for  one  of  the  nine-point  models  con¬ 
sidered  in  the  investigation.  Figures.  5  and  6  show  the 
effect  of  the  random  error  seed  on  a  twelve-point  model. 

All  computer  experiments  which  incorporated  error  used  a 
+5  percent  random  and  a  +5  percent  bias  on  the  real  and 
Imaginary  displacement  mobility  data. 

Figures  7,  8,  9,  10,  11  and  12  present  the  reidentified 
acceleration  mobility,  both  real  and  imaginary,  for  typical 
five-,  nine-,  and  twelve-point  models  respectively.  The 
models  varied  in  that  different  spanwise  masses  were  con¬ 
sidered.  Some  of  the  models  employed  in  the  study  are  given 
in  Table  VII  showing  the  various  points  of  interest  for  each 
model.  For  each  model,  the  computer  experiments  were 
executed  using  the  same  random  number  seed  and  the  afore¬ 
mentioned  errors  were  incorporated.  As  evidenced  by  the 
figures,  the  various  models  provided  acceptable  reidentifica 
tion  of  the  twenty-point  specimen  simulated  experimental 
displacement  mobility  data. 
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(POSITIVE) 
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Figure  1.  Effect  of  Error  on  Five-Point  Model 
Identification  of  Real  Acceleration 


FROM  IDENTIFIED 
MODEL"- _ 


3 

H 


H 

«  m 

OS 

o 

9 

«  <#> 

W  (0 

S  in 

« 

o  u  w 

U 

W 

2  >*  >• 

% 

2  £ 

Q 

E-«  O 

U 

H 

M  Q 

w 

n  oi 

2  2 

(0 

OH* 

(0  5 

w 

UJ 

o  <n  m 

w 

w  <#> 

< 

o\  o\ 

<  in. 

u 

fN  tN  fN 

U  + 

fN 

I 

O 


c 

o 

•H 

+> 

ns 


(aoHoa-ai/o)  ahvnidvwi  ' NoiivnaaaDDY 


- - — - — — iinmrtnr  mawiMiiifiiinri 


31 


Figure  2  .  Effect  of  Error  on  Five-Point  Model 
Identification  of  Imaginary  Acceler 
Response?  Driving  Point  at  Hub. 
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Figure  3.  Effect  of  Error  on  Nine-Point  Model 
Identification  of  Real  Acceleration 
Response;  Driving  Point  at  Hub. 


100 

Figure  5 .  Effect  of  Error  on  Twelve-Point  Model 
Identification  of  Real  Acceleration 
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Figure  9.  Effect  of  Model  on  Nine-Point  Model 
Identification  of  Real  Acceleration 
Response;  Driving  Point  at  Hub. 


Figure  10  -  Continued 
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Figure  11  -  Continued  . 
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Figure  12.  Effect  of  Model  on  Twelve-Point  Model 

Identification  of  Imaginary  Acceleration 
Response  ;  Driving  Point  at  Hub. 


CONCLUSIONS 


1.  Single-point  excitation  of  a  structure  yields  the 
necessary  mobility  data  to  satisfactorily  determine 
the  mass,  stiffness  and  damping  characteristics  for  a 
mathematical  model  having  less  degrees  of  freedom  than 
the  linear  elastic  structure  it  represents. 

2.  The  method  does  not  require  an  intuitive  mathematical 
model  and  uses  only  a  minimum  amount  of  impedance- type 
test  data. 

3.  The  eigenvector  or  mode  shape  associated  with  each 
natural  frequency  is  also  determined  in  the  analysis. 

4.  Computer  experiments  using  simulated  test  data  indicate 
the  method  is  insensitive  to  the  level  of  measurement 
error  inherent  in  the  state  of  the  measurement  art. 

5.  A  fully  populated  mass  matrix  should  be  assumed  for  an 
accurate  analytical  model  of  a  real  structure. 


i 


i 
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APPENDIX 

COMPUTER  PROGRAM  DESCRIPTION 


A  digital  computer  program  was  designed  for  computer  experi 
ment  to  investigate  the  proper  physical  interpretation  of 
identified  parameters  for  use  in  helicopter  engineering. 

The  program  was  written  for  the  IBM  360/40  operating  system 
using  FORTRAN  IV  language.  A  flow  chart  indicating  the 
program  logical  procedure  is  shown  in  Figure  13.  A 
description  of  the  input  cards  and  a  program  source  listing 
are  included  in  this  appendix. 


SUBROUTINE 


FORM  MODE 
SHAPE  MATRIX 
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PSEUDO 
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READ 
INITIAL 
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Flow  Chart  of  Computer  Program. 
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Figure  13  -  Continued. 
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INPUT  STRUCTURAL  DYNAMICS  PROGRAM  STIDN 
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SINGLE  POINT  FORCING 


C 

c 

c 


c 

c 

c 

c 

c 

c 


c 


STRUCTURAL  OVNANICS 


MNP 

MNP 


MNP 

INTEGER  HEAD! 201 ,HT(7I  MNP 

DIMENSION  SI20, 211 , VI (20,211  .3SUI2J),  MNP 

AYR  I  20(21  ),PHI(20,2n,PHtA(2P,21),UMF(20l,G(?Ot,SH20,21),  MNP 

dSNI20.21l.SII20.2i) ,751(20,211,  MNP 

CARS  I  201 , ANSI  201 ,A0S(2O) .0NNI20) , YRTI 2D , 100 > , YITI 20,  100)  MNP 

DIMENSION  PHINI20.21) , Y I Nl 20 ,21 ) , PH  I  Ml  20,21)  MNP 

DIMENSION  2SRI 20,21 1 ,OMFSI 2n  I , JMFKl 20 ) ,TR(lO0,20),TItl00,2/ll  MNP 

DIMENSION  OMNCI 20) , AKSRI 20)  ,  A HSR 1 20)  .ClNNS (2O),A0SR(20l  MNP 

DIMENSION  HZI100),INOXI20  I ,KEcP( 23)  ,  PHI T I  20,21 )  MNP 

DIMENSION  DPRIIOO,  20), OPIIIOO, 20), 2RI20, 21), 21  1  20,21)  MNP 

LOGICAL  TORF  MNP 

DATA  HT/*EXEC*,'T  DA'.'TA  S • , • MUL • , • ATED* , •  TES*.*T  •/  MNP 

DISPLACEMENT  M031LITY  DATA  ON  TAPE  13  1184)  MNP 

NJ-NUMBER  OF  GENERALIZED  COORDINATE)  MNP 

NUNUER  OF  OEGREES  OF  FREEDOM  MNP 

NQ- NUMBER  OF  MOOES  MNP 

NP*NUMBER  OF  FORCING  FREQUENCIES  MNP 

NK-FORCE  INPUT  STA  MNP 

REMIND  14  MNP 

READ  11,140)  IPl,IP2,NR0M,NN,NJ,NK,i TMS.NFF  MNP 

REAO  (1,140)  I  KEEP! I ) » I ■ l »N J  )  MNP 

130  READ  11,120)  ATOL ,  PTOL  ,PCTR,PCTilK,PCTI,PCTBI*IZ,IA  MNP 

I X* 1 2*2+1  MNP 

READ  (1,140  )  NPM1  MNP 

REMIND  13  Mf  P 

READ  (131  NCOL.HT, HEAD.NF.NO, I  iUI I )  ,1-1 ,NF  I  MM> 

DO  110  L-l.NF  1MNP 

UJ  READ  (13)  HZ(LI,((YRT(I,LI,YIT(I,L)I,1*1,ND  )  1MNP 

123  FORMAT  (4F10.4.2UO  t  MNP 

MRITE  13,130)  PCTR.PCTBk.PCTI.PCTBI.U  MNP 

liD  FORMAT  I • 1  * , T10 , 'MAX  RAND  ERROR  ON  ML.  '  a’F6. 3, lnX,*BIAS  ERROR  MNP 

A  ON  REAL  • • F6«  3 , *  OF  ELEME ITS* /T 10, *  ON  IMAGINAR  MNP 

dYa'F6.3,21X,'CN  I  MAGI  NARY"  •  F6.J, ISA, 'SEED*1'  I5///I  MNP 

14}  FORMAT  (8110)  "NP 

ISO  FORMAT  I8F10.4)  MNP 

MRITE  13,170)  IKEEPII ),I-1,NJ  )  MNP 

MRITE  (3,160)  KEEPINK)  MNP 

16J  FORMAT  I /////•  FORCE  INPUT  IS  AT  STA  *13/7/1  MNP 

173  FORMAT  I///// T10, • STAT IONS  USED  • // 1 120 , IP  15) I  MNP 

PTOL«PT0L/10O.  MNP 

00  480  MMal.NPHI  1MNP 

READ  11,140)  NQ.NP  1MNP 

READ  (1,150)  IOMFI  I I.I-l.NP)  1NNP 

DO  180  I « 1 , NP  2 MNP 

UMFRt  II «ONF(  I  1*6.283185  2 MNP 

180  ONFSI II-OMFRI !)*OMFR(I I  2MNP 

REAO  11,140)  I  INOX  I  11,1*1 «NP  )  IMNP 

READ  INITIAL  S  MATRIX  (ROMMISEI  IMNP 

DO  190  I *1 , NQ  2MNP 

19 j  REAO  (1,150)  I  SII  ,J) ,J-1 ,NP)  2MNP 

CALL  RE 01  I VR,Yt,NP,NJ,KEEP.INJA  , YRT, Yl T  )  IMNP 

1  TC  *  1  I  IMP 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 
3  7 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


56 


oo  uoou  ooooo  ou 


wRITfc  (  3.2001  ( CMF ( It . 1  *  1 , NP) 

IMNP 

56 

20J 

FORMAT  (////////// T50,  'PORCIN)  FRLJJtNCI ES  '//(10F12 

.4)/////)  IMNP 

57 

WRITE  (3.210) 

I  MNP 

58 

21) 

FORMAT  ( •  1  • , T50 , *R£  At  MORI L I T Y  MATRIX'//) 

l  MNP 

59 

CALL  MOUT?  (VR.NJ.NP  1 

IMNP 

60 

WRITE  (3,2201 

1MNP 

61 

22  D 

FORMAT  (*1',T50,  'IMAGINARY  MOIUITY  MATRIX'//) 

IMNP 

62 

CALL  MOUT 2  (YI.NJ.NP  I 

IMNP 

63 

IFtPCTR.NE.O.OR.PCTBR.NF.o.nR.PCTI.NE.O.OR.PCTBI.NE.O 

)  CALL  ERR  Nil  IMNP 

64 

i 

4  ( YR.VI  .PCTR.PCTBR.PCTI .PCTBl .NJ.NP, IX  1 

IMNP 

65 

WRITE  (3,2301 

IMNP 

66 

2)0 

FORMAT  ('l'»T50»' MOBIL  1  TV  MATRlCta  WITH  ERROR  REAL. 

IMAGINARY')  IMNP 

67 

CALL  MOUT 2  (VR.NJ.NP  1 

IMNP 

68 

CALL  MOUT?  (VI.NJ.NP  1 

IMNP 

69 

IMNP 

70 

NORMALIZE  IMAGINARY  MOBILITY 

IMNP 

71 

IMNP 

72 

IMNP 

73 

IMNP 

74 

ITERATE  FOR  MODE  SHAPE  ANO  S  MATRIX 

IMNP 

75 

24) 

CALL  PSEUDO  (S.NQ.NP.SM  1 

IMNP 

76 

WRITE  I3.2SO)  ITC 

IMNP 

77 

25) 

FORMAT  (•  S  ITERATION* • 141 

IMNP 

78 

CALL  MOUT 2  (S.NO.NP  1 

IMNP 

79 

IMNP 

80 

IMNP 

81 

CALL  MMPV  (  VI  .SM.NJ.NP.Na. PHI  I 

IMNP 

82 

IMNP 

83 

IMNP 

84 

2a) 

FORMAT  (///•  PHI  MATRIX'/) 

IMNP 

85 

IMNP 

86 

IMNP 

87 

IMNP 

88 

NORMALIZE  PHI  MATRIX 

IMNP 

89 

CALL  ANORM  (  PHI ,PHI M.NJ, NO  ) 

IMNP 

90 

CALL  PSEUDO  I  PHI ,N J, NO, PHI  A  ) 

IMNP 

91 

IMNP 

3? 

IMNP 

S3 

IMNP 

94 

IMNP 

95 

IMNP 

96 

27) 

CALL  MMPV  (  PHIA.VI  ,NQ,NJ,NP,Sl  1 

IMNP 

97 

IMNP 

98 

IMNP 

99 

CALL  TRAN  1  SI.SM.NQ.NP  1 

IMNP 

100 

CALL  ANORM  (SM.ST.NP.NO  1 

IMNP 

101 

CALL  TRAN  (ST, SI. NP, NO  I 

IMNP 

102 

CHECK  CONVERGENCE  OF  S  MATRIX 

IMNP 

10  3 

DO  300  1*1, NQ 

?MNP 

1"4 

DO  300  J-l.NP 

3MNP 

105 

DEL*  SI  (I.J)-S(I.J) 

3MNP 

106 

IF  (  ABS  (DEL )  -  ATOL  )  30", 30", 28) 

3MNP 

107 

2e  j 

IF  ($11, Jl)  2S0.31O.2QO 

3MNP 

108 

29J 

IF  ( ABS(OtL/S( 1 ,J) S-PTOL)  3O),J)J,J10 

3  MNP 

109 

JO) 

CONTINUE 

3  MNP 

110 

\ 


\ 


57 


CO  TO  360 

1MNP 

111 

313 

IF  1  ITC-I TMS  I  320,320,34" 

I  MNP 

112 

32  J 

ITC-ITC+l 

IMNP 

113 

00  33C  J«1,NP 

2MNP 

114 

00  330  |al,NO 

3MNP 

115 

33J 

3(1, Jt-  SIII.JI 

3MNP 

116 

CO  TO  240 

IMNP 

117 

3*3 

WRITE  13,350) 

IMNP 

lie 

353 

FORMAT  IT10, 'MAXIMUM  NUMBER  OF  i  MATRIX  ITERATIONS  EXCEEDED, 

IMNP 

119 

«JOb  TERMINATED' 1 

IMNP 

120 

GO  TO  8 TO 

IMNP 

121 

363 

MR1TE  (3,260) 

IMNP 

122 

CALL  M0UT2  (  PMIM.NJ.NO  ) 

IMNP 

123 

WRITE  (3,370  ) 

IMNP 

124 

373 

FORMAT  (•  CONVERGEO  S  MATRIX'/) 

IMNP 

125 

CALL  M0UT2  (SI.NO,NP  ) 

IMNP 

126 

c 

CALCULATE  MODAL  MOBILITY 

IMNP 

127 

c 

SM-Y*  REAL  SI-Y*  IMAG 

IMNP 

128 

CALL  PSEUDO  ( PHIM.N J,NQ, PHIN  I 

IMNP 

129 

CALL  MMPV  IPHIN.TR, NO, NJ,NP,  SN  1 

IMNP 

130 

CALL  MMPY  (PHIN,YI,NO,NJ,NP,  SI  ) 

IMNP 

131 

WRITE  13,3801 

IMNP 

132 

38  J 

FORMAT  I'l'.TlO, 'MODAL  MOBILITIES,  REAL,  IMAGINARY'//) 

IMNP 

133 

CALL  M0UT2  (  SM.NO.NP  1 

IMNP 

134 

CALL  M0UT2  I  SI  ,NO,NP  1 

IMNP 

135 

c 

IMNP 

136 

c 

IMNP 

137 

c 

CALCULATE  MOOAL  IMPEDANCE 

IMNP 

138 

c 

IMNP 

139 

00  390  l-l.NQ 

2MNP 

140 

WRITE  (3,150  )  PH I MINK,  I  > 

2MNP 

141 

DO  390  J-l.NP 

3MNP 

142 

C0N-PH1M1NK.I  ) / ( SI ( 1 ,  J I •  SI(I,J>»  34(1, Jl*  SHI  I , J ) ) 

3MNP 

143 

ZSRII.JI-  S'  (  1 , J)*CON 

3MNP 

144 

393 

ZSKI.JI—  SI  1 1  ,J)*CON 

3MNP 

145 

WRITE  13,400) 

IMNP 

146 

403 

FORMAT  (•!',  T10,  'MODAL  IMPEDANCE  REAL  ,  1  MAG  INARV //  J 

IMNP 

147 

CALL  M0UT2  I2SR  ,NQ,NP  ) 

IMNP 

148 

CALL  M0UT2  (ZSI.NO.NP  I 

IMNP 

149 

c 

IMNP 

150 

c 

IMNP 

151 

c 

LEAST  SQUARES  ANALYSIS  ON  HOUAL  IMPcOANCE  AS  FUNCTION 

IMNP 

152 

c 

OF  FORCING  FREQUENCY  SQUARED 

IMNP 

153 

c 

IMNP 

154 

c 

IMNP 

155 

NL-NP/NQ 

IMNP 

156 

ANL-NL 

IMNP 

157 

NLC-NL 

IMNP 

158 

M-l 

IMNP 

159 

DO  420  K*1,NQ 

2MNP 

160 

SUM  -0. 

2MNP 

161 

SUMA-R. 

2  MNP 

162 

SUMB-O. 

2MNP 

163 

SUMC-O. 

7MNP 

164 

00  410  I  »K J , NLC 

3MNP 

165 

58 


SUM  -OMFSIII*SUM 

3MNP 

166 

SUMA-ZSRIK.I I ♦SUM* 

3MNP 

167 

SUMB-ONFSII )*OMFSI 1 >*SUMB 

3MNP 

16R 

[ 

+uj 

SUMO  »OMFSIll*ZSR(  K,  1  USUMC 

3MNP 

160 

UET»ANL*SUMB-SUM*SUM 

2MNP 

170 

XA-I SUMA6SUM8-SUMC*SUMI/0ET 

2MNP 

171 

i 

XB*  1  ANL  *SUMC-SUMA*SUM 1 /OET 

2MNP 

172 

I 

KJ-NLC*  1 

2MNP 

173 

\ 

■ 

NIC-NL*IK*1I 

2MNP 

174 

} 

OMNCIKI -S0RT1ABS  IXA/XB  II 

2MNP 

175 

AKSKlKI --XB*OMNC(  K 1 *OMNC ( K  t 

2MNP 

176 

AMSHIK) — XB 

2  MNP 

177 

'  • 

UMNSIK) »OMNClKI*OMNC(KI 

2  MNP 

178 

+ij 

CONTINUE 

2MNP 

170 

L«1 

1MNP 

180 

DO  430  I » 1  ,N0 

2MNP 

181 

AUSH  I  II -(QMFSIO/OMNSI  I  1-1.  I*  SI  I  i , L 1 *AK SR (I  I / 

SMII.LI 

2MNP 

182 

ONNC ( II -OMNCI  1 1  /6. 2831  B 

2MNP 

183 

430 

L-2* 1 ♦! 

2  MNP 

184 

C 

1MNP 

185 

c 

IMNP 

186 

IF  (  MM.NE.l  1  GO  TO  450 

1MNP 

187 

SUM-0. 

IMNP 

188 

: 

00  440  1-1, Nl 

2MNP 

180 

44  3 

SUM-ZSII  I  « I  I ♦SUM 

2MNP 

i«>n 

G(  1  1 -SUM/I  AXSRI  1  1 ♦ ANL 1 

IMNP 

101 

OMNI  1 l-OMNCI 1 1 

IMNP 

102 

J 

AOSIll-AOSRIll 

IMNP 

103 

ANSIll-AMSRIll 

IMNP 

104 

AKSI  ll-AKSRUI 

IMNP 

105 

WRITE  1141  1  PH  III,  MM  I.I-l.NJ  1 

IMNP 

106 

GO  TO  480 

IMNP 

107 

(i 

00  460  I-l.NJ 

2MNP 

108 

463 

PHITII.MM  1  -  PHI  11.21 

2MNP 

100 

SUM-0. 

IMNP 

200 

> 

NI-NL+1 

IMNP 

201 

NZ-2  *Nt 

IMNP 

202 

00  4 TO  I-NI.NZ 

2MNP 

2"3 

473 

SUM-ZSI1  2, II ♦SUM 

2MNP 

204 

01  MMI -SUM/I AK  SR  1  2  l*ANLI 

IMNP 

205 

OMNI  MMI-0MNCI2I 

IMNP 

206 

AOS  I  MMI-A0SRI2I 

IMNP 

26  7 

AMSI  MMI-AMSRI2I 

IMNP 

708 

AKSI  MMI-AKSRI2I 

IMNP 

200 

WRITE  1141  IPHITII.MM  t,I-l,NJ  1 

IMNP 

210 

4b  3 

WRITE  1  3.540  1  MM.OMNIMMI .AMSIMNI  .AKSIMMI 

.AOSIMMI 

IMNP 

211 

WRITE  (3.4001  t  OMNI  1  1 , 1-1  .NPrII  1 

MNP 

212 

WRITE  (  3,500  1  (  AKSI II ,1-l.NPHil 

MNP 

213 

- 

WRITE  (3.510  |  (  AMSI 1 I.I-1.NPHII 

MNP 

214 

4»J 

FORMAT  I/////TK, 'CALCULATED  NATURAL  FREOUENCI 

ES,  CrCLES/SECV 

MNP 

215 

4 

1  IP 10£ 1 3*41 1 

MNP 

216 

533 

FORMAT  I/////T10, 'CALCULATED  GENERALIZED 

STIFF NESS* /I1P10E13. 211 

MNP 

717 

. 

513 

FORMAT  I/////110,  •CALCULATED  Gc.«lRAL  1 ZEU 

MASS' 

/( 1P10EI 3.21  1 

MNP 

718 

REWIND  14 

MNP 

210 

i 

DO  520  J-l.NPHI 

IMNP 

220 

, 

59 

i 

\ 

| 

\ 

} 

UMNSI Jl -OMNI JI*OMN( J) 

1  MNP 

221 

GSQI JI-GI J)«G(J> 

1  MNP 

222 

5/J 

HEAD  (14)  ( PH  II  I  >  J  )  1 1  ■  l  >  N J  1 

IMNP 

223 

NQ-NPH1 

MNP 

224 

CALL  ANURM  (PHI .PHI M.NJ,  NQ  ) 

MNP 

225 

WRITE  ( 3.530  ) 

MNP 

226 

533 

FORMAT  ('l'.T50. 'NORMAL  MOOES'//  ) 

MNP 

227 

CALL  MOUT2  (PHIM.NJ.NO  ) 

MNP 

228 

CALL  PSEUDO  1  PHIM.NJ.NO.  PHIA  1 

MNP 

229 

c 

MNP 

230 

c 

MNP 

231 

c 

IDENTIFICATION  OF  MASS. STIFFNESS  ANJ  DAMPING  MATRICES 

MNP 

232 

CALL  TRAN  (PHIA.PHIN.NO.NJ  1 

MNP 

233 

543 

FORMAT  (///•  MOOAL  PARAMETERS  MJOL',14//' 

NATURAL 

FREQUENC  MNP 

234 

A  V-'  F  14.  3,  •  HERTZ'//'  GENERALIZED  MASS  -'F14.3 

,'  SLUGS' 

//•  MNP 

235 

OOENERAL IZEO  ST1 FF- ' FI 4 . 2  ,  •  LB/IN'//'  GENERALI  ZED  DAMP 

■  •F14.2.  MNP 

236 

O'  Lb-SEC/IN' /////) 

MNP 

237 

c 

SM-INVERSE  OF  NASS 

MNP 

238 

c 

ST-INFLLENCE  COEFFICIENT 

MNP 

239 

c 

SI-INVERSE  OF  DAMPING 

MNP 

240 

DO  560  J-l.NJ 

1  MNP 

241 

00  560  K-l.NQ 

2  MNP 

242 

SUM  1-0. 

2MNP 

243 

SUMM-O. 

2  MNP 

244 

SUMO-0. 

2MNP 

245 

00  550  I-I.NQ 

3MNP 

246 

ACON-PHIMIK.I l«PHIMIJ,  1) 

3  MNP 

247 

SUM  I -ACON/AKS ( I ) *  SUM I 

3MNP 

248 

SUMM-ACON/AMSI I )*SUMM 

3MNP 

249 

550 

SUMO -AC  ON / 1 AK SI II*GII) >*SUMD 

3MNP 

250 

STIK.JI-SUMI 

2  MNP 

251 

SM(K.J) -SUMM 

2MNP 

252 

560 

SI  I  K. J) -SUMO 

2  MNP 

253 

CALL  INVRS  (SM.NJ.ZSR  ) 

MNP 

254 

WRITE  (3.570) 

MNP 

255 

573 

FORMAT  I 'l',T50, 'IDENTIFIED  MASS  MATRIX'//) 

MNP 

256 

CALL  M0UT2  1 Z  SR  .NJ.NJ  ) 

MNP 

257 

WRITE  (3.580  ) 

MNP 

258 

5b3 

FORMAT  I 'l',T50, 'IDENTIFIED  I NFLUcNCE  COEFFICIENT 

MATRIX'// 1  MNP 

259 

CALL  M0UT2  (ST, NJ.NJ  ) 

MNP 

260 

CALL  INVRS  (ST.NJ.7SR  ) 

MNP 

261 

WRITE  13,550) 

MNP 

262 

553 

FORMAT  I'1',T50,'I0ENTIFIE0  STIFFNESS  MATRIX'//) 

MNP 

263 

CALL  M0UT2  I  ZSR. NJ.NJ  1 

MNP 

264 

WRITE  ( 3.(001 

MNP 

265 

630 

FORMAT  1 c 1 • , T50, ' I OENT IF IED  DAMPING  MATKtX'//l 

MNP 

266 

CALL  INVRS  I  SI ,NJ. ZSR  ) 

MNP 

267 

CALL  MUUT2  (ZSR  .NJ.NJ  1 

MNP 

268 

SUM-0. 

MNP 

269 

DO  610  I-l.NQ 

IMNP 

270 

WRITE  I  3.620)  I, GUI 

IMNP 

271 

613 

SUM- SUM *G 1 1 ) 

IMNP 

272 

GS-SUM/NQ 

MNP 

273 

620 

FORMAT  (18, F22, 4) 

MNP 

274 

WRITE  (3,630)  GS 

MNP 

279 

60 


tffifwww. 


6  iU 

FORHAT  (//•  AVG  STRUCTURAL  DAMP  1 NG» • E8.  61 

MNP 

276 

IF  (NFF.EQ.O  1  GO  TO  650 

MNP 

277 

66  j 

HEAU  11.150  1  (HZOl.l-l.NFF  1 

MNP 

278 

NF-NFF 

MNP 

279 

GO  TO  660 

MNP 

280 

650 

IF  (NF.E0.0)  GO  TO  870 

MNP 

281 

660 

TORF-NROH.GT.  0.  ANO.  NROH.  LE.NQ 

MNP 

282 

00  750  L-l.NF 

1MNP 

283 

CON-HZI LI*HZ(L) 

IMNP 

286 

CALL  HOBPHI  1  G.GSQ.CON. AMS ,OMN&,  YR, VI » PHI M,NQ, NJ 

1 

1  MNP 

285 

6/3 

IF  1 1 P 11  680.680,730 

IMNP 

286 

6b  3 

IFIIP2.NE.0I  CALL  MATAMP  f H2< U . VR , Ti  .NOI 

1  MNP 

287 

1FIIP2.NE.0I  GO  TO  700 

IMNP 

288 

HRITE  13.6901  HZ(L» 

IMNP 

289 

693 

FORMAT  1 •  l 'T60, 'REAL  MOBILITY,  IMAGINARY  MOBILITY 

FREQ  -'F10.2, 

IMNP 

290 

< 

k  •  HERT  Z*  // 1 

IMNP 

291 

GO  TO  720 

IMNP 

292 

703 

HRITE  (3,7101  HZIL) 

IMNP 

293 

713 

FORNATI 'l'T60, 'ACCELERATION  AMPL1TUJE  IN  G"S,  PHASE  IN  OEG.  FREQ1MNP 

296 

A  ••  F 10.  2*  •  HERTZ'//I 

IMNP 

295 

720 

CALL  MOUT2  (YR.NO.NQI 

IMNP 

296 

CALL  MOUT2  (YI.NO.NQI 

IMNP 

297 

GO  TO  750 

IMNP 

298 

730 

00  760  I-l.NQ 

2MNP 

299 

OPRIL.I  l-YRII.I) 

2  MNP 

300 

OPIIL.I  1  ■VIII, II 

2MNP 

301 

IF  I .NOT .TORFI  GO  TO  760 

2MNP 

302 

IRIL.II-VRINROH.II 

2MNP 

303 

Til L, 11 ■YIINROM.lt 

2MNP 

306 

763 

CONTINUE 

2  MNP 

305 

753 

CONTINUE 

IMNP 

306 

IF! IP  1)  870,870.760 

MNP 

307 

760 

IFIIP2.NE.il  GO  TO  780 

MNP 

308 

CALL  AMP  IHZ.OPR.DPI.NF.NOI 

309 

IFtTORFI  CALL  AMP  IHZ, TR.TI , NF.NJI 

•'MNP 

310 

HRITE  13,7701 

/ 

MNP 

311 

773 

FORMAT  l'l*T60, 'DRIVING  POINT  RESPONSE,  AMP  IN  G1 

"S  ANO  PH  AS  f-VN 

MNP 

312 

ADEGREES'//) 

MNP 

313 

GO  TO  810 

MNP 

316 

7s3 

HRITE  13.7901 

MNP 

315 

703 

FORMAT  (• l'TAO,  'DRIVING  POINT  NUBILITY.  REAL  ANO 

IMAGINARY'//) 

MNP 

316 

IF  (  IA.NE.0  >  HRITE  (3,8001 

MNP 

317 

800 

FORMAT  IT60, 'ACCELERATION  MOBILITY'//! 

MNP 

318 

813 

CALL  VOUT  (HZ ,0PR,NF,NQ,0. 1  A  1 

MNP 

319 

HRITE  (3.8201 

MNP 

320 

823 

FORMAT  t'l'//l 

MNP 

321 

CALL  VOUT  (HZ.OPI,NF,NO,IP2,I  A  ) 

MNP 

322 

IF(  .NOT. TORFI  GO  TO  870 

MNP 

323 

IF  ( IP2.NE.il  GO  TO  860 

MNP 

326 

HRITE  (3,8301  NROH 

MNP 

325 

830 

FORMAT  ('l'T30,  'TRANSFER  RESPONSE,  ROM  '15,'  AMP 

IN  G"S  ANO  PHAS 

MNP 

326 

A 

E  IN  OEG'//l 

MNP 

327 

GO  TO  860 

MNP 

328 

860 

HRITE  (3,8501  NROH 

MNP 

329 

853 

FORMAT  t'l'T30,  'TRANSFER  MOBILITY,  XOH  *15,'  REAL 

ANO  IMAG* // 1 

MNP 

33C 

61 


IF  I1A.NE.0  1  WRITE  13,8001 

MNP 

331 

860 

CALL  TOUT  4H2 ,TR  ,NF  , NQ,0, 1 A  1 

MNP 

332 

WRITE  13,820) 

MNP 

333 

CALL  TOUT  (HlfTI ,NF,NQ,  IP2.IA  1 

MNP 

334 

a  tj 

CONTINUE 

MNP 

KEWINO  13 

MNP 

336 

CALL  EXIT 

MNP 

33T 

END 

MNP 

338 

I 


62 


m m 


SUBROUTINE  TRAN  1  A,R,  NR.NC  1 

TRN 

1 

c 

B-TKANSPOSE  OF  MATRIX  A 

TRN 

2 

c 

A-UNOISTURBED  MATRIX 

TRN 

1 

DIMENSION  A(20,21» , B<20,21> 

TRN 

A 

00  100  IM.NR 

1  TRN 

5 

DO  100  <laI«NC 

2TRN 

6 

130  BIJ. 11-All, J) 

2  TRN 

7 

RETURN 

TRN 

8 

END 

TRN 

9 

o  r> 


SUBROUTINE  INVRS  IB.N.AI 

INV 

l 

A  -  INVERSE  OF  B  8  UNO  1 jTJRdEO 

INV 

2 

INV 

3 

OIMENSION  At 20, 21 I , 0120, 21 1 .  IRJAI  2i 1  .1C0L 1 21 ) . Bl 20, 21 1 

INV 

4 

00  100  1*1, N 

1  INV 

5 

00  100  J*1,N 

2  INV 

6 

103 

Al l,JI*B( I .  J  > 

2  INV 

7 

N*N»1 

INV 

e 

00  110  1*1, N 

1  INV 

9 

1R0MIII-I 

1  INV 

10 

110 

ICOLI II ■! 

1  INV 

11 

00  26C  K*1,N 

1  INV 

12 

AMAX-  AIK,K) 

l  INV 

13 

00  130  1*K,N 

2  INV 

14 

00  130  J*K,N 

3  INV 

13 

IFIABSI  A (I.J)I-ABSI AMAX 11130 ,123, 120 

3  INV 

16 

uo 

AMAX*  A 1 1 ,  J I 

3  INV 

17 

1C- 1 

3  INV 

18 

JC-J 

3  INV 

19 

133 

CONTINUE 

3  INV 

20 

K1-1C0LIKI 

1  INV 

21 

1 COL  IK) ■ ICOLI 1C  1 

1  INV 

22 

icon  ici*»i 

1  INV 

23 

K 1*1  ROM  IK) 

l  INV 

24 

IROMl K) ■IROMl JCI 

1  INV 

25 

IROMIX  )*KI 

l  INV 

26 

IF! AMAX 1  160,140,160 

1  INV 

27 

143 

MRITE  13,130) 

l  INV 

28 

130 

FORMAT! '  SOLUTION  OF  EXISTING  MATRIX  NOT  POSSIBLE' I 

1  INV 

29 

GO  TO  330 

1  INV 

30 

160 

DO  170  J-l.N 

2  INV 

31 

E*A 1 K, J 1 

2  INV 

32 

AIK, JI-AI IC , J 1 

2INV 

33 

170 

AIIC.JI-E 

2  INV 

34 

00  180  1*1, N 

2  INV 

35 

E-AII.K) 

2  INV 

36 

A4 1 , KI-AI 1 , JC ) 

2  INV 

37 

UO 

All • JC) *E 

2  INV 

38 

00  210  1*1, N 

2  INV 

39 

IF1I-KI  200,190,200 

2INV 

40 

103 

All ,MI* 1, 

2  INV 

41 

GO  TO  210 

2  INV 

42 

203 

All , M)*0, 

2  INV 

43 

213 

CONTINUE 

2  INV 

44 

FVT*AIK,K) 

1  INV 

45 

00  220  J*1,M 

2  INV 

46 

223 

AIK,J)*A(K,J)/PVT 

2  INV 

47 

00  230  1*1, N 

2  INV 

48 

IF  I I-KI 230,250,230 

2  INV 

49 

230 

AMULT*A 1 1 ,KI 

2INV 

50 

00  240  J*1,M 

7  INV 

31 

260 

AII,JI*AII,J) -AMUL  T*A 1 K , Jl 

3  INV 

52 

230 

CONTINUE 

2INV 

53 

DO  26C  1*1,  N 

2  INV 

54 

263 

All, K) ■All, M) 

2  INV 

55 

64 


CPBMMMMMumjmimi 


00  290  1-1, N 

1INV 

56 

00  270  L-l.N 

2INV 

57 

1FI IRON! 1 l-L 12 70,280,270 

2 INV 

58 

273 

CONTINUE 

2 INV 

39 

2B3 

00  290  J-1,N 

2  INV 

60 

293 

OIL,  JI-AI 1, Jl 

2  INV 

61 

00  32C  J*1,N 

t  INV 

62 

00  300  L-1,N 

2  INV 

63 

1FIIC0LUI-U  300,310, 309 

2  INV 

6A 

303 

CONTINUE 

2  INV 

65 

313 

00  320  I-l.N 

2INV 

66 

323 

A< 1 ,L)*0( l,J» 

2  INV 

67 

333 

RETURN 

INV 

68 

END 

INV 

69 

SUBROUTINE  MHPV  i A . B.N1 « N2.N3 ,CI 

HPV 

c 

HPV 

c 

c  •  a  *  a 

HPV 

c 

A  INI  X  M2 >  B  (N2  X  M3 1 

C  INI  X  N3> 

HPV 

c 

HPV 

REAL  At  20*211 *B( 20 *21 1 *CI 20*211 

HPV 

00  100  l-l.Nl 

IHPV 

00  100  J«1*N3 

2HPV 

CII.J1-0. 

2HPV 

00  100  R-1.N2 

3HPV 

1 

100  CU.J)-C(I.J>*A(I*K>*B(K.JI 

3HPV 

1 

RETURN 

HPV 

1 

END 

HPV 

1 

66 


SUBROUTINE  H0UT2  (A,M,N) 

MOT 

1 

HEAL  A( 20,100) 

MOT 

2 

IO*MINO (N, ln) 

MOT 

3 

UNITE  <  3> lOOt  II. 1-1, ID) 

MOT 

4 

133 

FORMAT  (/T5, 10112) 

MOT 

5 

MITE  13,1001 

MOT 

6 

00  110  l>l,N 

1M0T 

7 

no 

WRITE  13,1201  I , I A ( I , J I ,  J»1 ,10) 

1M0T 

a 

U3 

FORMAT  C5,5X,1P10E12.4) 

NOT 

9 

IF  I ID-NI  130,170,170 

NOT 

10 

uo 

I0-MINOIN.2OI 

MJT 

11 

WRITE  13,1001  (I,I>U,IO  1 

MOT 

12 

WRITE  13,100) 

MOT 

13 

00  140  1*1, M 

1M0T 

14 

143 

WRITE  13,120)  l,IAII,JI * J»1 1,10  ) 

INOT 

15 

IFIIO-NI  150, 1 7C, 1 70 

NOT 

16 

153 

WRITE  (3,100)  (1,1 *21 , N  ) 

MOT 

17 

WRITE  (3,100  1 

MOT 

18 

00  160  (-1.M 

INOT 

19 

163 

WRITE  (3,120  )  1,1  A(I ,J),J*2l,N  ) 

IMOT 

20 

170 

RETURN 

MOT 

21 

END 

MOT 

22 

SUBROUTINE  ANORM  ( PHI f PHINtNR • 4C  ) 

NRM 

DIMENSION  PHI  120*21 1 • PHI N( 2ft, 211 

NRM 

00  120  !>1,NC 

1NRM 

AMAX«PH!(  1,11 

INRM 

00  100  J»2*NR 

2NRM 

IFIABSI  ANAXI.lE.ABSIPHIt  J,imA<MX»PHI(J'M 

2NRM 

103 

CONTINUE 

2NRM 

00  110  j-l.NR 

2NRM 

113 

PHINI J, 1 1  «PHI ( J ,1 1 /AMAX 

2NRM 

U3 

CONTINUE 

INRM 

10 

RETURN 

NRM 

11 

ENO 

NRM 

12 

68 


SUBROUTINE  EHRNU  1  A  ,8  ,  PC  TR,  PC  TjR.Pi.  T  1  ,PC  TBI  ,  NJ,NP,IX) 

ERR 

c 

ERR 

c 

*  BIAS  ERROA, 

ERR 

c 

PC  TB  (RATIO)  ON  ANPLI TJJfc,  AND  A  UNIPORN  RANOOM  ERROR 

ERR 

c 

HAVING  A  ♦/-  MAXIMUM  Op  PCI  (RATIO)  ON  AMPLITUOE. 

ERR 

c 

ERR 

c 

ERR 

c 

ERR 

c 

USES  RANOU 

ERR 

c 

ERR 

10 

DIMENSION  A ( 20 «21 )  «8(  2r‘«  21 1 

ERR 

11 

IMPCTRI  110,100,110 

ERR 

12 

uj 

If (PCT8RI  110,130, 1 10 

ERR 

13 

113 

00  120  I-l.NJ 

1ERR 

1 A 

00  120  J-1,NP 

2ERR 

15 

CALL  RANOU  (IX.IV.VfLl 

2ERR 

16 

IX*  iV 

2FRR 

IT 

E* 1 <0+2 ,OPPCTR* ( VFL-0.5I*PCTBR 

2  ERR 

18 

AII,JI>AII,JI«E 

2  ERR 

15 

CALL  RANOU  (IX,IV,VFLI 

2ERR 

20 

SX-IV 

2  ERR 

21 

E«1.0»2.0*PCTI  •( VFL-0.5)*PCTB( 

2ERR 

22 

U3 

b ( 1 , J 1 *B( I , J ) *E 

2ERR 

23 

133 

RETURN 

ERR 

2A 

ENO 

ERR 

25 

69 


SUBROUTINE  RANOU  (IX,IV,YFll 

RAN 

1 

c 

THIS  SUBROUTINE  IS  FROM  SSP  VERS.  II 

RAN 

2 

1 V«IXA65539 

RAN 

3 

IF  1  IV)  100,110,1  0 

RAN 

A 

mo 

1  V«  1  V*2 1A7A8  “l*» 

RAN 

<i 

110 

VFL* IV 

RAN 

6 

VFL»VFL9.Ai, 

RAN 

7 

RETURN 

RAN 

8 

ENO 

RAN 

9 

SUBROUTINE  REOI  1 YR , VI  ,NP, NJ.KEEP ,UOX  .VRT.VIT  1 

RAN 

10 

c 

RAN 

11 

c 

REDUCES  DISPLACEMENT  MOBILITY  DATA  TO  MATRIX  OF  NJ  SPCCIMEN 

RAN 

12 

c 

COORDINATES  AND  FORCING  FREQUENCIES  V*NJ«NP 

RAN 

1) 

DIMENSION  VRI20.21  ),VI<20,21  » ,KEePI20> , INOXI 20  1 

RAN 

1 A 

DIMENSION  VRTI20.100I  ,V!T(2O,U0t 

RAN 

IS 

DO  120  I»1 ,NP 

IRAN 

16 

DO  120  J-l.NJ 

2RAN 

17 

VRIJ.I  )>VRTIKE£P(JI,'NOXU  II 

2RAN 

IB 

uo 

VliJtl l»YITIXEEPI  Jl > INOXI 1 1) 

2RAN 

19 

RETURN 

RAN 

20 

ENO 

RAN 

21 

70 


SUBKOUT INE  YOUT  <OMH, A,N INC .NO.^AHP , I A  1 

YOT 

1 

YOT 

2 

IF  1A  NOT  -  0  USE  ACCELERATION  ABILITY 

TOT 

3 

YOT 

4 

REAL  0MH(100),Ail0Q,2O| 

YOT 

5 

IF  I  IA  1  100.120,100 

YOT 

6 

133 

CON-  6.201185*6.283185 

YOT 

7 

00  110  l-l.NINC 

1  YOT 

8 

QM-OMHI 1 1 *OMH 1 1 1 *C  ON 

1  YOT 

9 

00  110  3-1, NO 

2V0T 

10 

113 

Ail  ,31— A 1 1 ,3I*CM 

2  YOT 

11 

Uu 

31-1 

YOT 

12 

10-NIN0IND, 101 

YOT 

13 

no 

IL-MIN0(NINC,45) 

YOT 

14 

11-1 

YOT 

15 

143 

WRITE  13,1501  (1,1-Jl.IOI 

YOT 

16 

153 

FORMAT  CT5,  •HERT2M6, 91121 

YOT 

17 

WRITE  13,1601 

YOT 

18 

163 

FORMAT  11X1 

YOT 

19 

1FINAMP 1  170,170,200 

YOT 

20 

173 

00  160  1-11, IL 

1  YOT 

21 

163 

WRITE  13,1901  OMNI  II, (All, 3) ,J-31, 1 Jl 

1  YOT 

22 

193 

FORMAT  <1X,F9.3,1P10E12.4) 

YOT 

23 

CO  TO  2  30 

YOT 

24 

233 

00  210  1-11, IL 

1Y0T 

25 

210 

WRITE  13, 2201  OMMII 1,1  At  1,31 ,J-Jl,l3» 

1V0T 

26 

223 

FORMAT  (1X,F9.3,10F12.2) 

YOT 

27 

230 

IF!  IL-NINC)  240,260,260 

YOT 

28 

263 

WRITE  13,2501 

YOT 

29 

253 

FORMAT  C  •  1  • // • 

YOT 

30 

11-46 

YOT 

31 

1L-NIIC 

YOT 

32 

GO  TO  140 

YOT 

33 

263 

IFIIO-NOI  270,280,200 

YOT 

34 

273 

Jl-U 

YOT 

35 

10-N0 

YOT 

36 

WRITE  13,220) 

YOT 

37 

GO  TO  130 

YOT 

38 

263 

RETURN 

YOT 

39 

ENO 

YOT 

40 

71 


SUBROUTINE  AMP  (OMH, A ,8 , NI NC, NKI 

AMP 

1 

c 

AMP 

2 

c 

CONVERTS  A  ♦  1*8  IN  DISPLACEMENT 

UNI  TS 

AMP 

3 

c 

TO  AMP  (IN  A  1  IN  G'i  AND  PHASE  (IN 

8  1  IN  OEG 

AMP 

6 

c 

EACH  ROW  IS  AT  A  MfcJUtNCV  OMHI II 

IN  HERTZ 

AMP 

5 

c 

AMP 

6 

c 

.0162  c. t'83ie5  /  386. 

AMP 

7 

DIME  Zl  -*  i  >jMH ( 100 >,A(  100 ,20  1,8(103.201 

AMP 

8 

c 

AMP 

9 

00  210  I-l.NINC 

1  AMP 

10 

ON-OMHI 11*0.01626 

(AMP 

11 

OMRaOMHIl  1*6.283185 

1  AMP 

12 

00  210  J-1'NR 

2AHP 

13 

R-A( l,JI 

2AMP 

16 

C*B( 1 « J 1 

2  AMP 

15 

A(ltJ)«SQRT(ft*AtC*CI*OK*OMR 

2  AMP 

16 

IF ( R  )  160.100,160 

2  AMP 

17 

1,3 

1FICI  110,120,130 

2AMP 

18 

113 

Bit , J)a  270. 

2  AMP 

19 

CO  TO  210 

2AMP 

20 

Uj 

B(I,JI-0 

2AMP 

21 

GO  TO  210 

2  AMP 

22 

133 

8(1 ,JI»90. 

2AMP 

23 

GO  TO  210 

2AMP 

26 

143 

P-ATAN(  ABS(C/RI  1*57.2458 

2AMP 

25 

IFIKI  150,130.180 

2  AMP 

26 

153 

IFICI  160,16n,l70 

2  AMP 

27 

163 

B(I,JI-180.*P 

2  AMP 

28 

GO  TO  210 

’’AMP 

29 

173 

BII.JI-180.-P 

2AMP 

30 

GO  TO  210 

2AMP 

31 

U3 

IFICI  140,190,200 

2  AMP 

32 

150 

8(1 , JI-360.-P 

2AMP 

33 

GO  TO  210 

2AMP 

36 

230 

BII.Jl-P 

2AMP 

35 

213 

CONTINUE 

2AMP 

36 

RETURN 

AMP 

37 

END 

AMP 

38 

72 


SUBROUTINE  CINV  |A,B,N,C,D> 

C1N 

1 

c 

CIN 

2 

DIMENSION  *120,211  , 8(20,21 1  , Cl  23,21 1 .01 20, 211, El?*,  ?tl 

CIN 

1 

c 

C*l'  «  INVtrtSt  OF  A»I*B  l-SORTI-n 

CIN 

A 

c 

CIN 

5 

c 

8  ASSUMED  NUN  SINGULAR 

CIN 

6 

c 

CIN 

7 

CALL  lNVRSie,N,CI 

CIN 

R 

CALL  MNPV(C«A,N,N«N,E) 

CIN 

9 

CALL  MNPT(A,E,N,N,N,CI 

CIN 

10 

00  100  !*1,N 

1C  IN 

11 

DO  100  J-1,N 

2C  IN 

12 

130 

CII.JI-CII.JI+BII.J) 

2CIN 

1) 

CALL  IN VRSIC, N,OI 

CIN 

1 A 

CALL  NNPVIE,D,N,N,N,CI 

CIN 

15 

DO  110  l«l,N 

1C  IN 

16 

DO  110  J«1,N 

?C  In 

17 

11J 

0<l,JI--0(I,JI 

2CIN 

18 

RETURN 

CIN 

19 

END 

CIN 

20 

I 


1 JBROUT INC  M08PHI  (  G, GSU.CON , A MS , iMNS , YR , V  1 , PHI* ,NO,N J  1 

MOB 

l 

c 

CALCULATES  VP  AND  VI  USING  W'lOAL  HUdlLlTV  AND  MOOE  SHAPE 

HOB 

2 

DIMENSION  G 1  201  ,GSQI20 1 ,  AMSI20I ,  VRI 20 .20  )  .  VI 1 20,201 , PHIMI  20,20), 

MOB 

3 

AVSRI20)  ,  VS  1 1 201  tOMNSI  26) 

MOB 

A 

00  100  1"  1,NQ 

1N0B 

5 

CONA-CON/OMNSII) 

1M0B 

6 

CONd- 1.  /  I  CON*  AM  SI  1  1*39.4  76413  i 

IMOB 

7 

CONC-CONA-1. 

1  MOB 

R 

CQNO-CONA  *CONB/ 1  CONC*CONC  »GSO  1 1 1 1 

IMOB 

9 

VSR  1 1 )a-CONC*CONO 

IMOB 

10 

103 

VSIIII— G(l)«CONO 

IMOB 

11 

00  120  J-l.NJ 

IMOB 

12 

00  120  K-l.NO 

2M0B 

13 

SUMR-O. 

2M0B 

1A 

SUM 1-0. 

2M0B 

15 

00  110  1-1, NO 

3M0B 

16 

ACON-PHIMIK,  1  IPPHIMI  J,  1  1 

IMOB 

17 

SUMR-VSRI  1 l*ACON»SUMR 

IMOB 

ia 

11  j 

SUM  1 -VS  II 1 )*ACON«-SUMI 

3M0B 

19 

VRIK.JI-SUMR 

2M0B 

20 

120 

Vl 1 K , Jl -SUMI 

2M0B 

21 

RETURN 

MOB 

22 

END 

MOB 

23 

74 


SUBROUTINE  PSEUDO  <A,NR,NC,CI 

PSU 

1 

c 

PSU 

2 

c 

C«-  PSEUDOINVERSE  OF  A  A  UN01STURBF0 

PSU 

3 

c 

A  IS  A  RECTANGULAR  MATRIX  JF  MAXIMAL  RANK  INR 

X  NCI 

PSU 

A 

c 

NR  .GT.  OR  .LT.  NC 

PSU 

5 

c 

PSU 

6 

c 

-1  -l 

PSU 

7 

c 

C  -  (A'A)  A*  OR  AMAA'I 

PSU 

8 

c 

PSU 

9 

c 

NR, NC  NAY  NOT  EXCEEO  25 

PSU 

10 

c 

PSU 

11 

REAL  A(20.21)«B(29,2L)<C(20,21I 

PSU 

12 

c 

li  •  A* 

PSU 

13 

UO  100  I-l.NR 

IPSU 

1 A 

00  100  J-l.NC 

2PSU 

15 

130 

Bl J,  II-AI  I,JI 

2PSU 

16 

I F  (  NR-NC ) 120,110,1 30 

PSU 

IT 

113 

CALL  INVRS  I A  ,  NR ,  C  1 

PSU 

18 

GO  TO  1  AO 

PSU 

19 

c 

NR  ,LE.  NC 

PSU 

20 

c 

C  ■  AA* 

PSU 

21 

1x3 

CALL  MMPY  ( A , B , NR , NC , NR , C ) 

PSU 

22 

c 

A  -  INV  UF  C 

PSU 

23 

CALL  INVRS  (C ,NR, A) 

PSU 

2A 

c 

C  -  PSEJUJINVEKSE 

OF  A  (NC  X 

NR  1 

PSU 

25 

CALL  MMPY  IB, A, NC, NR, NR, Cl 

PSU 

26 

GO  TO  1  AO 

PSU 

27 

c 

NC  .LT.  NK 

PSU 

2ft 

c 

C  *  A*  A 

PSU 

29 

133 

CALL  MMPY  (B,A,NC,NR,NC,CI 

PSU 

30 

c 

A  -  INV  JF  C 

PSU 

31 

CALL  INVRS  IC.NC.AI 

PSU 

32 

c 

C  «  PSEUJJl  INVERSE 

OF  A  INC  X 

NR) 

PSU 

33 

CALL  MMPV  ( A  ,B, NC , NC , NR , C 1 

PSU 

3A 

c 

RESTORE  A 

PSU 

35 

UJ 

00  150  I-l.NR 

1°SU 

76 

DO  150  J-l.NC 

2PSU 

77 

1!>J 

All  .JI-BIJ.II 

2PSU 

3ft 

RETURN 

PSU 

79 

ENO 

PSU 

A  r> 

75 


